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NOMENCLATURE 


a  Isentropic  speed  of  sound 

h  Static  enthalpy 

H  Ratio  of  local  coordinate  system  curvature  to  body  radius 

of  curvature,  1  +  y/R 

Total  (stagnation)  enthalpy,  h  +  l/2q2 

K  Curvature  of  body 

M  Mach  number 

p  Pressure 

P  Logarithm  of  pressure 

q  Total  velocity,  / u2+v2+w2 

r  Radial  coordinate  in  cylindrical  coordinates 

Cross-section  radius 

R  Body  radius  of  curvature,  or  gas  constant 

s  Entropy 

T  Temperature 

u  Tangential  velocity  component  in  body-normal  coordinates 

v  Normal  velocity  component  in  body-normal  coordinates 

V  Freestream  velocity 

oo 

w  Circumferential  velocity  component  in  body-normal 

coordinates 

x  Wetted  length,  measured  along  body  surface  from 

stagnation  point 

y  Distance  from  body  surface,  measured  along  body-normal 

i  Axial  coordinate  in  cylindrical  coordinates 
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NOMENCLATURE  (Cont'd.) 

a  Angle  of  attack 

a*  Stretching  parameter 

6  Sideslip  angle 

y  Isentropic  exponent 

e  Convergence  criterion  for  iterative  calculations 

n  Transformed  normal  coordinate 

0  Circumferential  coordinate  in  (^n,0)  coordinates 

0j3  Local  body  angle 

X  Characteristic  slope 

C  Wetted  length  in  U,n,@)  coordinates 

p  Density 

a  Normal  streamline  slope,  v/u 

as  Local  shock  slope,  tan'^rj/az 

t  Circumferential  streamline  slope,  w/u 

<f>  Circumferential  coordinate 

$  Circumferential  coordinate  (nose  geometry) 

{  ^  Condition  at  body 

(  )EB  Equivalent  body  quantity 

(  )s  Condition  at  shock 

(  )0  Stagnation  point  condition 

(  1  Freestream  condition 

'00 

(  )*  Sonic  point  condition 
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SECTION  1 


INTRODUCTION 

This  report  describes  work  performed  under  Task  3. 2. 1.9,  "Asym¬ 
metric  Flow  Field  Modeling,"  Subtask  3. 2. 1.9. 2,  "3DSAP  Modifications,"  of 
the  Performance  Technology  Program  (PTP-S  II).  This  effort  involved  the 
extension  of  the  Three-Dimensional  Shock  and  Pressure  (3DSAP)  approximate 
Inviscld  flow  field  technique  to  treat  flared  multi-conic  frustum  con¬ 
figurations. 

The  3DSAP  approximate  flow  field  technique  was  originally  de¬ 
veloped  to  provide  rapid  yet  accurate  aerodynamic  predictions  for  ballistic 
reentry  vehicles  with  asymmetric  ablated  nosetips  (References  1  and  2). 

In  3DSAP  the  nosetip  flow  field  solution  is  developed  by  assigning  surface 
pressures  from  a  correlation  and  iterating  on  bow  shock  position  until 
mass  continuity  is  satisfied,  assuming  the  functional  forms  of  the  profiles 
of  the  flow  properties  across  the  shock  layer.  The  frustum  flow  field 
solution  is  obtained  from  the  finite  difference  solution  of  the  three-dimen¬ 
sional  inviscid  flow  equations,  using  initial  data  from  the  nosetip  flow 
field  solution.  The  solution  of  the  3-D  inviscid  equations  on  the  frustum 
,»  simplified  by  assuming  the  functional  form  of  the  cross-flow  distribution 
and  limiting  the  analysis  to  small  angles  of  attack  and  axi;ymmetric  frusta. 
(.,s  originally  developed,  the  3DSAP  code  was  restricted  to  conic  or  expan¬ 
sion  biconic  frusta.) 

Section  2  of  this  report  describes  an  optional  nosetip  pressure 
correlation  that  has  been  added  to  3DSAP  for  use  on  indented  nosetip 
geometries.  This  optional  correlation  was  developed  in  an  earlier  task 
of  the  PTP-S  II  program  for  use  on  indented  nosetip  geometries  in  a 
dispersion  code  (Reference  3). 


The  modification  of  the  3DSAP  code  to  allow  radial  grid  stretching 
in  the  supersonic  finite  difference  calculations  is  described  in  Section  3. 
The  addition  of  this  option  to  the  code  permits  the  clustering  of  mesh 
points  near  the  wall  to  resolve  the  flow  gradients  that  may  arise  there, 
while  retaining  the  formal  second  order  accuracy  of  the  finite  difference 
procedure  used. 

In  Section  4  the  extension  of  the  geometric  capabilities  of  the 
3DSAP  code  to  multi-conic  axisymmetric  afterbodies  is  described.  Included 
in  this  new  capability  is  an  Improved  procedure  for  treating  the  disconti¬ 
nuities  of  both  the  flow  and  the  coordinate  system  at  expansion  corners. 

Finally,  in  Section  5,  the  addition  of  the  capability  of  treating 
flared  afterbody  sections  to  the  3DSAP  code  is  described.  In  this  analysis 
it  is  assumed  that  there  is  an  attached  embedded  shock  on  the  flare  and 
that  the  flow  behind  the  shock  remains  supersonic.  A  discrete  shock-fitting 
algorithm  has  been  developed  to  accurately  treat  the  embedded  shock. 

The  appendices  to  this  report  provide  a  modified  user's  guide 
to  the  3DSAP  code. 
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SECTION  2 


MODIFIED  NOSETIP  SURFACE 
PRESSURE  CORRELATION 

As  described  in  Reference  1,  the  3DSAP  approximate  flow  field 
procedure  assigns  nosetlp  surface  pressures  from  correlations,  assuming 
that  each  meridional  plane  on  the  nosetlp  Is  treated  as  an  equivalent 
axlsymmetrlc  body.  In  the  original  3DSAP  analysis,  the  approximate  loca¬ 
tion  of  the  sonic  point  Is  determined  from  a  correlation  relating  the  sonic 
point  body  angle  to  the  Inverse  bluntness  of  the  nosetlp,  as  developed  In 
Reference  4.  In  the  subsonic  region  upstream  of  the  sonic  point,  the  sur¬ 
face  pressure  is  assigned  from  Love's  "Newtonian  Deficiency  Method,"5  as 
modified  in  Reference  4.  Downstream  of  the  sonic  point,  the  surface  pres¬ 
sure  is  assigned  from  a  curve  fit  developed  In  Reference  4  of  a  matched 
Newtonian  plus  Prandtl-Meyer  expansion. 

In  Reference  3,  as  part  of  the  PTP-S  II  program,  this  surface 
pressure  correlation  procedure  was  modified  to  provide  more  accurate  re¬ 
sults  for  indented  nosetlp  geometries  for  application  to  a  ballistic  reentry 
vehicle  (3RV)  dispersion  code.  This  Improved  correlation  has  been  In¬ 
corporated  as  an  option  In  the  3DSAP  code,  and  is  described  below. 

Using  the  nomenclature  of  the  original  3DSA15  equivalent  body 
analysis  (Reference  1),  the  sonic  point  body  angle  Is  defined  as 

8*  =  50°  -  aEp  (2.1) 

where  the  equivalent  angle  of  attack  is  defined  as 

aEQ  =  cos  0  +  B  sin  $  .  (2.2) 


7. 


The  wetted  length  from  the  stagnation  point  to  the  sonic  point  Is  denoted 
by  s*,  and  the  local  cylindrical  radius  Is  r*.  The  nosetip  pressure  cor¬ 
relation  is  applied  in  two  parts:  one  upstream  of  the  sonic  point  and  the 
other  downstream.  Upstream  of  the  sonic  point,  the  correlation  consists 
of  a  melding  of  correlations  for  spheres,  cones,  and  flat-faced  disks.  In¬ 
cluding  curvature  effects  for  Indented  shapes. 

The  spherical  pressure  correlation,  a  modification  of  the  form 
given  by  Kyriss  and  Neff",  is 

P, s  ■  1  -  1*08  (1  -  )  cos2  0bEQ  (2.3) 

where  p  *  p/p0,  Po  is  the  stagnation  pressure,  and  the  equivalent  body 
angle  is 


0bEQ  a  eb  +  aEQ  •  (2.4 

The  conical  pressure  correlation,  which  has  been  developed  in 
Reference  3  to  include  effects  of  concave  geometries,  is  given  by 


pc  *  min  {  1.0,  (1.15  sin20bEq  +  <$)  x 


(2.5) 


where  the  term  6  is  0  for  convex  geometries.  For  concave  geometries,  6 
accounts  for  curvature  effects  by  taking  the  form 


<5  *  cos  0bEQ  (r  2-  rj2) 


(2.6) 


where  ri 
and  K  is 


is  the  cylindrical  radius  at  the  start  of  the  concave  section, 

an  average  curvature,  defined  as 
s 


+  2K  ] 


(2.7) 
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(Here,  Is  the  wetted  length  at  the  start  of  the  concave  section.)  K 
is  the  local  curvature  at  a  given  point. 

The  spherical  and  conical  pressure  correlations  are  combined 
with  a  correlation  for  flat  disks  to  yield 

p  »  T?J  (  ps  -  ( ps  -  P  J  f FD 1 

♦  (1  -  V)  [pc  -  (pc  -  p*)  f  FD]  (2.8) 

where,  from  References  5  and  6, 

fFQ  =  B(s/s*)2  +  ( 1-B)  e'X  (2.9) 

with 

B  * - ^ - 

16  (1  - 

A  =  5  /in  (s*/s) 
p*  *  0.532 

The  weighting  factor,  IT,  is  defined  as 

IT  *  min  (1.0,  -  (tt/2  -  e^/s)  (2.10) 

where  Kq  is  the  curvature  at  the  stagnation  point. 

Downstream  of  the  sonic  point,  the  same  basic  approach  is  used, 
except  that  the  spherical  pressure  correlation  is  replaced  with  the  curve 
fit  of  a  matched  Newtonian  plus  Prandtl-Meyer  expansion  developed  by  Kyriss 
and  Neff4: 
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Ps  *  (1  -  jr  )[  1  *  1*^6143  sir,2  ($EQ  +  A6) 

+  0.51143  sin3  ( +  A0)]  +  ^  (2.11) 

where 

^EQ  =  f  '  0bEQ 
A0  =  0*  -  50° 

In  addition,  the  weighting  factor  R  is  redefined  to  be 


with  R*  =  r*/cos  0*.  Finally,  the  conical  pressure  term,  given  by  Equation 
(2.5),  is  constrained  to  not  exceed  the  term  p$  given  by  Equation  (2.11), 
and  the  resulting  pressure  is  given  by 

p  =  R3  p$  +  (1  -  R3)  p^  .  (2.13) 

This  optional  nosetip  pressure  correlation  is  invoked  auto¬ 
matically  by  3DSAP  for  non-snherical  nosetips,  or  can  be  selected  by  the 
user  by  setting  the  input  parameter  IPRESS  to  2.  Comparisons  of  the 
nosetip  surface  pressures  predicted  with  this  correlation  to  experimental 
data  for  indented  nosetips  may  be  found  in  Reference  3. 


10. 


SECTION  3 


NORMAL  GRID  STRETCHING  FOR 
SUPERSONIC  CALCULATIONS 

In  the  original  30SAP  analysis  of  Reference  1,  supersonic  calcu¬ 
lations  were  performed  in  a  computational  transformed  coordinate  system 
(E,n,8),  where  the  transformed  normal  coordinate  n  was  defined  as 

n  =  y/ys  (3.1) 

where  y  is  the  coordinate  normal  to  the  body  surface  and  y  =  ys  ( x ,d?)  de¬ 
fines  the  bow  shock  surface.  With  this  transformation  n  varies  between  0 
at  the  body  and  1  at  the  bow  shock,  and  computational  grid  points  equally 
spaced  in  n  are  also  equally  spaced  in  y  along  a  given  line  x  =  constant, 

<J)  =  constant. 

To  allow  for  radial  grid  stretching  the  coordinate  transformation 
defined  by  Equation  (3.1)  has  been  generalized  to  the  form 

n  =  f(y/ys)  •  (3.2) 

By  a  suitable  choice  of  the  function  f,  it  is  then  possible  to  cluster  mesh 
points  near  the  wall  which  are  equally  spaced  in  the  transformed  coordinate 
r\  (i.e.,  the  mesh  points  will  be  unequally  spaced  in  y).  Such  a  transforma¬ 
tion  will  permit  the  accurate  resolution  of  strong  flow  gradients  near  the 
wall  (e.g.,  entropy  layers)  while  retaining  the  formal  second  order  accuracy 
of  the  integration  procedure  used,  for  which  mesh  points  equally  spaced  in 
the  computational  coordinates  are  required. 

The  form  of  the  stretching  function  incorporated  into  3DSAP  is 
that  developed  by  Moretti  in  Reference  7.  In  this  formulation  the  stretched 
coordinate  n  is  defined  implicitly  from 

y/ys  =  1  +  tanh  (  a*(n-l)]/tan  h  a*  (3.3) 
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where  a*  is  the  stretching  parameter.  The  effect  of  this  formulation  on 
the  distribution  of  grid  points  is  shown  in  Figure  3.1;  note  that  as  the 
stretching  parameter  approaches  zero  the  effect  of  the  transformation  on 
the  grid  point  distribution  diminishes. 

With  the  coordinate  stretching  defined  by  Equation  (3.3),  the 
transformation  of  derivatives  in  the  physical  (x,y,40  coordinates  to  the 
computational  (C»n»6)  coordinates  becomes 


_3_  _  X 
3x  "  3^ 


+ 


(3.4) 


_3_ 

ay 


(3.5) 


x  =  x  + 

3<t>  30 


(3.5) 


where 


f  -  D 

y  ‘  a*[0My-ys)!] 

fx  =  fy 

%  ■  fy  ^  Wys 

D  =  y$/tan  h  a* 

Note  that  as  these  forms  reduce  to 
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* 


fx *sx/ys 
U  *  •"  J'sj./J's 

which  are  the  fonts  given  by  Equations  (2.2.9)  -  (2.2.11)  in  Reference  1, 
in  which  there  is  no  coordinate  stretching. 

With  coordinate  stretching  the  three-dimensional  inviscid  equa¬ 
tions  of  motion  become,  using  the  nomenclature  of  Reference  1, 

CONTINUITY 

(1  '  &)  pc  +  (A  -  Jr  fx)  Pn  +  BP0  +  yH  [fy 
52 

+  ~  <Vn  +  Te  +  t2  eb)  (3.7) 

6 1 

+  —  (sin  0b+o  cos  05)]  =  0 

MOMENTUM 

°£  +  %  +  Boe  *  l-°pc  +  (Hfy  -  afx>  p^i 

ltt>  (3-8) 

*  - tB  (COS  9(3  -  o  Sin  0b)  =  0 


+  Atd  +  Bt6  +  +  (Hf/r  -  Tfx)  Pn  +  y.  p0] 

(3.9) 

+  8  (sin  +  0  cos  0^ )  -  +  tzB  sin  9^  =  0 
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ENERGY 


s£  +  Asn  +  Bs6  s  0 

where 


(3.10) 


+  oHf..  +  Bf . 

y  4> 


B  =  62  tH/r 


At  body  points,  the  form  of  the  characteristic  compatibility  equa 
tion  that  is  solved  with  coordinate  stretching  takes  the  form 

=  -  XPn  -  [BPq  +  Y{onfy  +  (62Te  +  sin  6b)/r 

fw  cos  9k  f w  *2  (3.11 

+  tB  (sin  9b.  JL_b)  -  x{n/(l  -  jjr  ) 

with 

X  =  -fy/  / u2/a2-l 

At  bow  shock  points  the  new  characteristic  compatibility  condi¬ 
tion  takes  the  form 


^  61  <p«  *  *pn>  -  i  <A  -  f*)(“5  +  %> 

+  f  fy  <'5  *  Xvn>  '  0 


(3.12) 
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where 


R2  =  -  [  Bu0  +  ^  -  tuB  sin  eb  ] 

R3  =  "  [  Bve  "  ^  "  tuB  cos  eb  1 

In  the  modified  version  of  3DSAP,  the  Equations  (3.7)  -  (3.12) 
are  solved  exactly  as  were  their  counterparts  without  coordinate  stretching 
in  the  original  3DSAP,  with  only  one  small  variation.  It  has  recently  been 
noted  that  in  the  bow  shock  point  solution  procedure,  where  the  bow  shock 
slope  is  iterated  until  Equation  (3.12)  is  satisfied  to  within  some  tolerance 
e,  that  the  equation  is  in  a  dimensional  form.  Thus,  to  insure  an  accurate 
solution  to  the  bow  shock  equation,  the  convergence  criterion  is  made  con¬ 
sistently  dimensional  by  requiring  the  residual  of  Equation  (3.12)  to  be 
less  than  e/(Ax)2,  where  Ax  is  the  step  size  for  the  supersonic  calculation 
and  e  is  the  convergence  criterion  input  to  the  code  via  TEST(34).  This 
modification  ensures  a  valid  shock  point  solution,  regardless  of  the  length 
scale  used  in  the  geometric  input  to  3DSAP. 
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SECTION  4 

TREATMENT  OF  MULTI-CONIC  AFTERBODIES 


In  the  original  version  of  3DSAP  it  was  assumed  that  the  after¬ 
body  was  axisymnetric  ^id  was  either  a  cone  or  a  bicone.  To  extend  the 
afterbody  geometry  capability,  the  modified  3DSAP  allows  as  many  as  five 
conical  segments,  each  of  which  must  be  axisymmetric.  (The  number  of  coni¬ 
cal  segments  allowed  in  the  new  version  of  3DSAP  may  readily  be  increased 
by  increasing  the  dimensions  of  the  appropriate  arrays  within  the  3DSAP 
code. ) 


With  the  new  version  of  3DSAP  the  frustum  geometry  is  always  in¬ 
put  in  the  same  format.  As  shown  in  Figure  4.1,  the  I—  conical  segment 
is  defined  by  specification  of  RF(I),  ZF(I),  DZF(I),  THF(I) ,  RF(I+1),  and 
ZF( 1+1 ) .  Assuming  that  RF(1)  and  ZF(1)  are  known,  each  conical  segment  I 
is  defined  through  specification  of  any  one  of  the  following  pairs  of  para¬ 
meters: 


ZF( 1+1) ,  RF( 1+1 ) 

ZFO+l),  THF(I) 

RF( 1+1) ,  DZF(I) 

RF( 1+1) ,  THF(I) 

DZF(I),  THF(I) . 

The  choice  of  which  pair  of  parameters  to  use  is  left  to  the  user,  and 
can  be  varied  from  one  conical  section  to  the  next.  More  specific  details 
on  the  geometry  definition  procedure  for  the  new  version  of  3DSAP  are  pro¬ 
vided  in  Appendix  A. 

In  addition  to  increasing  the  number  of  frustum  geometry  segments 
allowed,  the  3DSAP  code  has  been  further  modified  to  improve  the  treatment 
of  expansion  corners  on  the  frustum.  Since  3DSAP  is  formulated  in  a  body- 
normal  coordinate  system,  the  coordinate  system  is  discontinuous  at  points 
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where  the  body  slope  Is  discontinuous.  As  illustrated  in  Figure  4.2,  the 
discontinuous  coordinate  system  actually  fails  to  include  a  portion  of 
the  flow  field  at  expansion  corners. 

To  illustrate  the  procedure  developed  to  treat  frustum  expansion 
discontinuities,  referring  to  Figure  4.2,  assume  the  flow  field  solution 
has  been  computed  to  the  end  of  the  first  conic  section  (cone  angle 
with  the  bow  shock  located  at  y  *  ys^,  and  with  the  flow  field  properties 
P,  o,  t,  s  known  along  the  body  normal  between  the  body  and  the  shock.  For 
the  forward  marching  supersonic  solution  to  be  computed  on  the  second  conic 
section  (cone  angle  602)  it  is  necessary  to  know  the  bow  shock  location 
y  *  ys2  measured  normal  vo  the  second  conic  section,  as  well  as  all  flow 
field  properties  along  this  new  normal. 

Because  of  the  coordinate  system  discontinuity,  no  flow  field 
information  is  available  along  the  new  body  normal.  However,  the  required 
information  can  readily  be  generated  through  the  use  of  several  approxi¬ 
mations. 


First,  it  is  assumed  that  the  bow  shocx  slope  as  measured  in  a 
continuous  (z,r)  cylindrical  coordinate  system  Is  constant  across  the  body- 
normal  coordinate  system  discontinuity*  Defining 


as  *  tan 


->  . 


9b,  +  tan^ysx 


(4.1) 


and  assuming  as  to  be  constant,  it  can  be  shewn  that 

cos  9bj  +  tan  os  sin  9b^ 

*2  *  ysl  cos  9(^2  +  tan  °s  cos  eb2 


(4.2) 


♦This  procedure  assumes  that  the  bow  shock  is  not  highly  curved  in  the 
region  of  the  coordinate  system  discontinuity. 
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from  which  it  follows  that 


rs2  »  rb  +  y$2 
ZS2  8  *b  *  ys2 

The  shock  slope  in  the  body-normal  coordinate  system  relative  to  the  second 
cone  can  then  be  expressed  as 

yS)<2  *  tan  (°s  *  ®b2)  .  (4.5) 

For  the  definition  of  flow  field  properties  along  the  new  normal, 
define  nj  =  f(y/ysj)  as  the  computational  coordinate  along  the  normal  to 
the  first  cone,  and  n 2  *  f(y/ys2)  as  the  computational  coordinate  along 
the  normal  to  the  second  cone.  Assuming  that  the  function  f  (which  may 
involve  stretching,  as  described  in  Section  3)  does  not  vary,  the  pressure 
and  entropy  are  assigned  as 

p( ^2 )  =  P(nl)  (4.6) 

s(n2)  3  s(nj)  .  (4.7) 

Note  that,  except  for  the  body  surface  point,  these  relations  would  hold 
exactly  if  the  first  conic  section  were  a  sharp  cone. 

More  care  is  required  in  the  modification  of  the  streamline  slopes 
a  and  x  due  to  the  discontinuity  of  the  body-normal  coordinate  system.  From 
purely  geometric  considerations,  it  can  be  shown  that  the  streamline  slope 
a2  on  the  second  cone  can  be  expressed  in  terms  of  the  slope  aj  on  the  first 
cone,  assuming  a  constant  streamline  slope  in  cylindrical  coordinates,  as 

Oj  +  tan  (9bj  -  9b2^  (4.8) 

2  1  -  tan  (6bx  -^b2) 


c°s  eb2  (4.3) 

sin  0b2  .  (4.4) 
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To  recompute  the  circumferential  streamline  slope  x  it  is  first 
assumed  that  the  total  velocity,  defined  as 

q  ■  / uj2  ( 1  +  Oj2  +  xj2 )  *  /ug*  (1  +  og2  +  T22)  (4.9) 

can  be  determined  from 


q(n2)  *  q(ni) 


(4.10) 


Since  the  cross-flow  velocity  component  w  is  not  affected  by  the  body-normal 
system  discontinuity,  it  follows  that 


T2 


u2  u2  Uj  u2 


(4.11) 


Combining  (4.9)  and  (4.11)  yields 


/q*  -  Tt~u7 
7  l  +  o22 


and  finally 


x2 


(4.12) 


(4.13) 


The  procedure  described  above  is  used  to  generate  the  necessary 
data  along  the  normal  to  the  second  cone  for  all  grid  points,  except  for 
the  body  point.  At  the  body  point  at  the  discontinuity,  the  pressure  on 
the  second  cone  is  redefined  from  the  series  expansion  the  Prandtl-Meyer 
relation  for  an  ideal  gas: 


YM  2 

7^r (Ae)  +  *i* 


(y+DMj4  -  4  (Mj2-  1) 


4(M12  -l)2 


(A6)  = 


(4.14) 
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where  A0  «  9k  -a.  TkJ. 

tions  to  approx  Invite*  the  pre^e'droVa^ss  thT  ^  T  *“ 

the  upstream  value  of  y.  th  ?dns<on  corner,  using 


condition  requires  th  Z  -  0 "j IIT"  ***'  ** 

quires  S,  ,  s„  rj «M"*ton  process  re- 

the  conservation  of  tota,  ^^“^"“vls'cidfl'o"^'  C°”PUted  ^ 


u2  *  ■''“ou  -  h2)  -  u,*  T;; 


(4. IS) 


the  body  cen^e^efi’ned^'^uaUoMTnT'"1''’  s'°pe  at 


SECTION  5 


CALCULATION  OF  FLARED  AFTERBODIES 

The  major  development  effort  in  the  modification  of  the  3DSAP 
code  has  been  the  addition  of  the  capability  to  treat  flared  afterbody 
sections.  The  calculation  procedure  for  flared  afterbody  sections  described 
in  this  section  was  developed  subject  to  the  following  assumptions: 

•  the  frustum  has  only  one  flare  (i.e.,  compression 
corner) ; 

•  the  flare  section  must  be  the  last  frustum  conic 
section; 

•  the  flare  section  flow  field  must  have  an  em¬ 
bedded  shock  of  finite  strength  that  is  attached 
at  all  points  of  the  compression  discontinuity; 
and 

•  the  flow  behind  the  embedded  shock  must  remain 
sufficiently  supersonic  so  that  the  f4nite- 
difference  forward  marching  integration  procedure 
for  supersonic  flows  remains  valid. 

With  these  assumptions,  a  discrete  shock-fitting  approach  was 
selected  as  the  solution  procedure  for  the  flared  afterbody  flow  field. 

In  this  approach  the  flare  shock  layer  is  divided  into  two  regions,  one 
between  the  body  and  the  embedded  shock  and  the  other  between  the  embedded 
shock  and  the  bow  shock.  The  embedded  shock,  which  serves  as  the  boundary 
between  the  two  regions,  is  explicitly  treated  as  a  discontinuity  across 
which  the  Rankine-Hugoniot  shock  jump  conditions  are  enforced. 

An  alternate  approach  to  the  calculation  of  the  flare  flow  field 
is  the  use  of  the  "conservation"  form  of  the  governing  equations,  which 
has  the  capability  of  automatically  "capturing"  the  embedded  shock  without 
requiring  any  special  procedures  to  calculate  the  position  and  strength 
of  the  embedded  shock.  However,  the  use  of  this  "shock-capturing"  technique 
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will  produce  oscillations  in  the  flow  field  in  the  vicinity  of  the  embedded 
shock;  these  oscillations  can  be  large  enough  to  cause  the  solution  to 
fail.  "Conservation"  form  solutions  frequently  use  numerical  damping  to 
control  these  oscillations  and  thus  obtain  a  solution,  but  if  too  much 
damping  is  used,  the  calculated  ir, viscid  flow  field  can  be  distorted. 

Because  of  this  consideration,  and  the  desire  to  have  an  algorithm 
that  will  work  with  minimal  intervention  required  of  the  user,  the  shock¬ 
fitting  approach  was  selected  for  this  effort  even  though  the  development 
of  a  shock-fitting  procedure  requires  greater  analytical  effort  than  does 
the  development  of  a  "shock-capturing"  scheme. 

The  shock- fitting  computational  scheme  has  been  incorporated  into 
a  new  subroutine,  FLARE,  which  is  called  from  the  subroutine  SHPR  when  a 
compressive  geometric  discontinuity  is  encountered  in  the  frustum  flow  field 
calculation.  (Thus,  when  using  3DSAP  on  configurations  that  do  not  have 
flares,  some  core  storage  can  be  saved  by  not  loading  the  subroutine  FLARE.) 

The  first  step  in  the  flared  body  calculation  is  the  initializa¬ 
tion  of  the  embedded  shock  at  the  compression  corner.  In  each  meridional 
plane  being  computed,  it  is  necessary  to  determine  the  slope  of  the  embedded 
shock  such  that  the  downstream  flow  is  tangent  to  the  body  surface.  Fol¬ 
lowing  the  procedure  outlined  in  Section  4  for  expansion  corners,  the 
computed  flow  field  solution  on  the  conic  section  upstream  of  the  flare  is 
rotated  to  the  body-normal  coordinate  system  for  the  flare  section,  but 
without  modifying  the  surface  pressure  or  imposing  the  kinematic  boundary 
condition  at  the  body  point.  Let  these  data  at  the  body  be  denoted  by 
Pl»  aj,  and  Sp  noting  that  aj  <  0. 

The  embedded  shock  initialization  problem  for  an  axisynmetric 
frustum  reduces  to  finding  the  slope  of  the  embedded  shock  such  that  the 
properties  downstream  of  the  shock  (denoted  by  P2,  T2,  and  S2)  satisfy 

the  flow  tangency  condition,  02  =  0. 


25. 


To  allow  for  real  gas  solutions,  an  iterative  process  is  used  to 
determine  the  initial  embedded  shock  slope  in  each  meridional  plane  being 
computed.  If  the  unit  normal  to  the  embedded  shock  is  denoted  by 


V 


N  ■  Ni  ex  +  N2  ey  (5.1) 

as  shown  in  Figure  5.1,  the  upstream  velocity  normal  to  the  shock  is 

ui  =  -Uj  Nj  -  u1  a1  N2  .  (5.2) 

Knowing  Qj,  p^,  and  sj,  the  downstream  properties  u2,  P2>  and  s2  can  be 
determined  from  the  Rankine-Hugoniot  conditions  from  a  call  to  the  3DSAP 
subroutine  SHOCK  for  both  ideal  and  real  gases. 

The  velocities  downstream  of  the  embedded  shock  can  be  determined 
from  the  relations 


u2  =  Uj  -  Nj  (Qj  -  U2)  (5.3) 

u2  cr2  s  Uj  oi  -  N2  (uj  -  u2)  .  (5.4) 

In  the  iteration,  the  value  of  N2  is  varied,  with 

Nx  =  -  /l-  N2z  ,  (5.5) 

until,  from  Equation  (5.4),  the  tangency  condition  a2  =  0  is  satisified. 
With  the  embedded  shock  surface  defined  as  y  =  ys2  ( x , ) ,  the  appropriate 
embedded  shock  slope  becomes 

yS2  ■  -N1/N2 

X 


(5.6) 


The  initialization  for  the  two  layer  flare  flow  field  calculation 
is  carried  out  as  shown  in  Figure  5.2.  With  the  bow  shock  surface  denoted 
by  y  =  ys1  (x,<j>),  Region  I  is  defined  as  yS2  <  y  <  ys^.  and  Region  II  as 
0  <  y  <  ys2-  Computational  transformations  of  the  normal  coordinate  are 
carried  out  separately  in  each  region  (and  may  Include  stretching),  result¬ 
ing  in 


=  f(S?)  (5'7) 

n2  =  f2(^)  (5‘8> 

where  the  form  of  the  functions  fj  and  is  that  given  in  Section  3  when 
stretching  is  used. 

The  initial  data  surface  is  established  at  a  small  but  finite 
distance  Ax  from  the  compression  corner  in  each  meridional  plane  as  shown 
in  Figure  5.2,  where  the  bow  and  embedded  shock  positions  are  given  by 


=  yn0  +  ysix  (4x) 

(5.9) 

ys2  =  ys2x 

(5.10) 

The  field  data  in  Region  I  is  assigned  directly  from  the  solution  upstream 
of  the  flare,  rotated  into  the  body-normal  coordinate  system  of  the  flare 
section.  This  region  will  initially  have  the  same  number  of  grid  points 
in  the  normal  direction  as  did  the  upstream  solution.  Region  II,  because 
of  its  small  size,  is  initialized  with  two  grid  points,  one  at  the  body 
and  one  at  the  embedded  shock,  and  both  points  have  the  properties  down¬ 
stream  of  the  embedded  shock  assigned. 
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The  computational  procedure  used  in  this  flared  body  subroutine 
is  exactly  that  used  in  the  subroutine  SHPR,  with  the  following  exceptions: 

1. )  The  "body"  point  in  Region  I  (which  is  the  upstream 

point  at  the  embedded  shock)  is  treated  as  a  field 
point,  using  one-sided  forward  differences  to  approxi¬ 
mate  n-derivati ves . 

2. )  Shock  point  calculations  (at  both  the  bow  and 

embedded  shocks)  are  performed  using  the  non¬ 
iterative  Kentzer-Moretti  predictor-corrector8*9 
algorithm. 

In  the  Kentzer-Moretti  approach  to  shock-point  calculations, 
differentiated  forms  of  the  Rankine-Hugoniot  conditions  are  incorporated 
into  the  appropriate  characteristic  compatibility  condition  to  yield  an 
equation  for  shock  curvature,  which  can  be  integrated  twice  to  obtain 
shock  slope  and  position. 

The  form  of  the  compatibility  condition  used  is  exactly  that 
given  by  Equation  (3.12),  with  the  additional  simplification  of  H  =  1 
allowed  by  consideration  of  only  conic  sections  on  the  afterbody. 

For  both  the  bow  and  embedded  shocks,  let  (  )^  denote  properties 

upstream  of  the  shock  and  (  ^  denote  downstream  properties.  It  is  then 

possible  to  differentiate  the  Rankine-Hugoniot  conditions  to  obtain 

p2^  =  D1  pl^  +  °2  Ql£  +  d3  sl£  (5.11) 

q2^  =  °4  Pl^  +  °5  al£  +  d6  sl^  »  (5*12) 
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) 


where 


Ui  =  "Ul  (Nj  +  Cj  N2  +  Tj  N3) 


Ni  =-ysx/y 


N2  =  1/v 
n3  = 

^  =  /  1  +  ySx2  +  ys^2/rs2 

It  can  be  noted  that  for  the  bow  shock,  where  Pj_  -  Jin  p^  and 
=  s^,  that  P^  =  =  0.  For  an  ideal  gas,  the  coefficients  Dj  may  be 

expressed  in  relatively  simple  form  as 

-  (Y-1)  a:2 
Dl  2YUJ2  -  (y-1)  ax2 

1  4yQi2 

°2  =  ^  ZyUj*  -  (y-1) 

-Ziij2  (y-1) 

°3  =  -  (y-1)  a^ 

2  (Y-1) 

°4  =  Y  (y+1)""  Uj 
°5  =  &  (Y'!  •  ) 
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For  a  real  gas  (equilibrium),  the  forms  of  the  coefficients  are  more 
complicated  and  are  not  presented  here. 


Other  terms  necessary  for  this  shock  point  equation  are  as 

fol lows: 


NU  =  Cl  y$xx  *  C2 

(5.13) 

"25  ‘  C3  ysxx  +  C4 

(5.14) 

N35  =  C5  ysxx  *  C6 

(5.15) 

“1;  *  c7  ysxx  *  c3 

(5.16) 

u25  ■  c9  ysxx  +  C10 

(5.17) 

v25  =  CU  ysxx  *  C12 

(5.18) 

coefficients  Ci  given  by 

C3  *  VNj 

C4  *  VN3  (ySx^/rs  -  yH  yH/rs>) 

C2  =  '*sx  C4 
C5  *  "ys4>  C3/rs 

C6  =  ^  Vrs  -  N2  -  *s*  ^Sx/rsM 
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C i  *  "U^  (Cj  +  Oj  Cg  +  Cg) 

Cg  *-u^  ( N1  +  Oj  N2  +  Tj  Ng)  -  Uj  (Ng  ai^  +  Ng  t^) 

“U1  ^C2  +  °1  C4  +  T1  C6^ 

Cg  ~  “(tig  ~  Uj)  Cj^  ~  (  Dg  “  1)  Cy 

C10  “  >Jlr  -«!  (D4  P1c  +  06  slc)  -  Nt  (Ds  -  1)  Cg 

-(Qg  ~  Uj)  Cg 

Cji  =  ~Ng  ( Dg  -  1 )  C?  -  (Qg  -  Uj )  Cg 

C12  *  °1  ul£  +  U1  alj;  ~  ^2  ^D4  Pl£  +  °5  +  C3  ^°5  ” 

“(Qg  -  0^)  C4 

With  substitution  of  the  above  expressions  into  the  characteristic 
compatibility  condition,  solving  for  the  shock  curvature  yields 

ysxx  =  (  0^  6i  (°i  Pl£  +  D3  s1^  +  u-  fx)  C10 

'  u  fy  C12  *  <RHS)1/'^  B1  D2  C7  (5.19) 

-  i  <fl  ~  f*>  C9  +  K  Cll] 

where  unsubscripted  quantities  are  evaluated  with  properties  downstream  of 
either  the  bow  or  embedded  shock,  as  appropriate. 
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Integrating  this  equation  twice  yields  both  the  shock  slope  (ysx) 
and  the  shock  position  (ys).  Knowing  the  shock  slope,  the  properties  down¬ 
stream  of  the  shock  can  be  determined  directly  using  the  SHOCK  subroutine. 

As  part  of  the  flared  afterbody  calculation  procedure,  logic 
has  been  included  to  automatically  respace  Regions  I  and  II  to  ensure  that 
sufficient  field  points  are  used  to  describe  the  region  between  the  body 
and  the  embedded  shock  as  the  shock  moves  away  from  the  body,  and  to  ensure 
that  the  mesh  spacing  between  the  embedded  and  bow  shocks  does  not  get  too 
small.  The  respacing  criteria  built  into  the  code  thus  ensure  adequate 
resolution,  without  unnecessarily  small  step  sizes  (which  will  adversely 
affect  the  efficiency  of  the  calculation). 

A  procedure  has  also  been  developed  to  treat,  in  an  approximate 
manner,  the  intersection  of  the  embedded  and  bow  shocks.  The  two  shocks 
are  assumed  to  intersect  when,  in  any  meridional  plane,  the  shocks  come 
within  some  arbitrary  percentage  of  tne  total  shock  layer  thickness  of 
each  other  (typically  5%).  When  this  occurs  in  any  one  plane,  the  shocks 
are  assumed  to  intersect  in  all  planes  at  points  determined  by  extending 
the  bow  and  embedded  shocks  with  constant  slopes  (ysx)  assumed.  At  the 
assumed  points  of  intersection,  the  exact  solutions  fo»  properties  down¬ 
stream  of  intersecting  shocks  of  the  same  family  are  obtained  (including 
the  slope  of  the  resultant  shock).  The  slopes  of  the  resultant  shocks  are 
then  extrapolated  back  to  the  last  station  computed,  a  new  "resultant  shock 
layer  thickness"  is  defined,  and  the  one  region  afterbody  solution  is  con¬ 
tinued  using  the  subroutine  SHPR. 

It  should  be  noted  that  in  many,  if  not  most,  of  the  flared  con¬ 
figurations  of  interest  that  the  bow  and  embedded  shocks  typically  inter¬ 
sect  downstream  of  the  aft  end  of  the  vehicle. 
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APPENDIX  A. _ USER'S  GUIDE  TO  3DSAP 


This  Appendix  describes  the  Input  variables  required  for  the 
operation  of  the  updated  3DSAP  code.  These  input  variables  are  discussed 
in  groupings  of  geometry  specification,  freestream  conditions,  output  con¬ 
trols,  and  special  options.  A  glossary  of  all  input  variables  is  pro¬ 
vided  in  Appendix  B. 

A.  GEOMETRY 


The  geometry  specifications  required  in  3DSAP  are  described  be¬ 
low.  Nosetip  geometries  may  be  arbitrary,  or  one  of  several  analytic 
options.  The  frustum  geometry  inputs  allow  for  the  definition  of  axi sym¬ 
metric  multi-conic  configurations. 

L  NOSE  GEOMETRY 


Arbitrary 


Arbitrary  nose  geometries  are  specified  by  defining  cross-sections 
perpendicular  to  the  vehicle  centerline,  as  shown  in  Figure  1.  Cross- 
sections  are  defined  in  terms  of  the  polar  coordinates  RBD(K,M)  and  PHIBD(M) 
at  the  axial  stations  ZBO(K),  relative  to  the  center  of  each  cross-section, 
defined  from  the  arrays  D1(K)  and  D2(K).  Note  that  the  cross-section  cen¬ 
ters  need  not  lie  along  the  vehicle  axis. 

Options  are  provided  for  specification  of  special  types  of  cross- 
sections,  as  discussed  in  Reference  2.  The  option  exercised  is  determined 
by  the  number  of  meridional  planes  of  geometry  input,  and  the  appropriate 
values  of  PHIBD  are  generated  automatically  as  follows: 
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No.  of  Planes  Input 
1 
2 
3 

4( ISA  =  2  only*) 


Geometry 
circular 
ell iptic 
bl-elliptlc 
quad-ell Iptic 


PHIBD 

0 

0,tt/2 

0,tt/2,it 

0,tt/2,  tt,  3tt/2 


If  more  than  four  planes  of  gemoetry  are  input  (three  planes,  for  ISA  s  1) 
the  cross-sections  are  taken  to  be  arbitrary  and  the  PHIBD  array  must  be 
specifi ed. 


The  surfaces  to  be  computed  in  the  nose  calculation,  defined  as 
4>  *  constant,  are  specified  by  the  user  with  the  PHI(M)  array.  Appropriate 
radii  at  each  cross-section  for  each  $  value  are  determined  analytically 
for  the  special  cross-section  options  listed  above;  for  the  case  of  general 
cross-sections  the  radii  are  determined  by  cubic  spline  Interpolation. 

If  the  user  does  not  specify  the  number  of  surfaces  to  be  computed, 
the  PHI  array  is  automatically  set  to  provide  for  surfaces  spaced  90°  apart; 
for  a  pitch  plane  of  geometric  symmetry  without  sideslip  (ISA  =  1),  the 
surfaces  to  be  computed  are  $  =  0,  tt/2,  and  tt.  Without  a  plane  of  symmetry 
(ISA  =  2),  the  surfaces  are  $  *  0,  tt/2,  tt  and  3tt/2.  ($  is  defined  such 
that  $  *  0  corresponds  to  the  lee  plane  for  positive  angle  of  attack,  and 

A 

4>  *  tt  to  the  wind  plane.) 

For  geometries  with  a  pitch  plane  of  symmetry  that  have  general 
cross-sections  and  are  to  be  run  at  sideslip  (requiring  ISA  =  2),  it  is 
possible  to  define  only  the  half-planes  of  the  arbitrary  cross-sections, 
imposing  symmetry  of  the  geometry  about  the  pitch  plane  by  specifying  the 
<i>  =  m  plane  as  MSYM. 


♦ISA  is  a  flag,  set  to  one  if  a  pitch  plane  of  gemoetric  symmetry  exists 
and  the  sideslip  is  zero;  otherwise,  ISA  is  set  to  two. 
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The  stagnation  point  on  arbitrary  nose  geometries  in  3DSAP  is 
defined  by  the  user-supplied  values  of  Dl(l)  and  D2(l);  RBD(l)  must  always 
be  zero.  This  point  will  be  common  to  all  surfaces  PHIBD(M)  and  PHI(M). 

In  addition,  the  last  point  of  the  nose  geometry,  which  must  coincide  with 
the  start  of  the  frustum,  must  have  D1  and  D2  both  zero,  and  all  values 
of  RBD  at  this  station  must  be  equal  to  RF(1)  (see  frustum  geometry  de¬ 
finition). 

Anal tyic 


Options  exist  within  the  3DSAP  for  the  analytic  definition  of 
spherical,  elliptic,  and  biconic  nosetips.  A  spherical  nosetip,  with 
radius  RN,  is  assumed  if  no  ZBD’s  are  input,  and  IB IC  =  0  and  IELL  =  0. 
The  axial  location  of  the  geometric  stagnation  point  is  assigned  to  be 
Z  =  DELZ. 


If  IELL  f  0,  an  elliptic  nosetip  is  assumed,  with  the  ratio  of 
major  to  minor  axes  defined  as  AB(a/b),  as  shown  in  Figure  A.l.  Setting 
IB  I C  f  0  allows  a  biconic  nosetip  to  be  defined  by  the  parameters  THN,  RC, 
and  RN,  as  shown  in  Figure  A.l.  Note  that  Z  =  0  is  defined  as  the  virtual 
apex  location  of  the  first  conical  afterbody  segment  for  both  the  elliptic 
and  biconic  nosetip  analytic  geometry  options. 

2.  FRUSTUM  GEOMETRY 


The  definition  of  the  frustum  geometry  in  this  modified  version 
of  3DSAP  requires  the  specification  of  values  for  the  arrays  ZF,  RF,  DZF, 
and  THF  (see  Figure  4.1).  The  frustum  is  assumed  to  be  an  axisymmetric 
multi-conic,  with  NCON  conic  sections  (maximum  of  5).  Values  for  RF(1) 
and  ZF(1)  must  always  be  provided  (except  for  a  spherical  nosetip),  and 
each  conical  segment  I  is  defined  through  specification  of  any  of  the 
following  pairs  of  parameters: 
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ZF(I+1),  RF( 1+1) 
ZF(  1+1) ,  THF(I) 
RF( 1+1) ,  OZF(I) 
RF(I+1),  THF(I) 
DZF(I),  THF(I) 


Cylindrical  segments  (THF(I)  =  0.)  cannot  be  defined  using  the  RF(I+1), 
THF(I)  pair. 

B.  FREESTREAM  PROPERTIES 


For  each  calculation  it  is  necessary  to  specify  the  freestream 
Mach  number  (AMINF)  or  the  freestream  velocity  (QINF),  as  well  as  the  static 
pressure  and  temperature  (PINF  and  TINF).  Values  of  PINF  and  TINF  will  be 
generated  automatically  from  the  1962  standard  atmosphere  tables  by  speci¬ 
fying  the  desired  altitude  (ALT)  and  the  flat  IATMP  =  62  for  the  1962  U.S. 
Standard  tables  or  IATMP  =  59  for  the  1959  ARDC  tables.  For  non-dimensional 
ideal  gas  calculations,  PNIF  and  TINF  are  set  to  1.0  (default  values). 

Care  must  be  exercised  when  specifying  QINF  for  non-dimensional  ideal  gas 
calculations,  since  then  QINF  =  Vw/ /p J p^",  as  discussed  below. 

The  orientation  of  the  vehicle  relative  to  the  freestream  velocity 
vector  is  defined  by  specification  of  the  angle  of  attack  (ALPHA)  and  the 
sideslip  angle  (BETA).  Note  that  positive  angle  of  attack  implies  nose 
up,  and  positive  sideslip  angle  implies  nose  left. 

Units 


IDEAL  GAS  ( IRG  =  0) 

The  standard  operating  procedure  for  ideal  gas  calculations  is 
to  use  non-dimensional  variables,  in  which  PINF  =  TINF  =  RIDEAL  =  1.  With 
this  normalization  the  variables  in  3DSAP  take  the  forms: 
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3DSAP 

P 

P 

T 

(u.v.w) 


ACTUAL 

P/P» 

P/P» 

T/T 

oo 

(u,v.w)/^75l 


Note  that  with  this  formulation  the  gas  law  becomes  p  *  pT.  The  isentropic 
exponent  is  specified  as  GIDEAL.  If  desired,  a  compressibility  factor 
(ZIDEAL)  may  also  be  specified. 


Dimensional  ideal  gas  calculations  are  also  possible,  requiring 
only  that  a  dimensionally  consistent  value  of  RIDEAL  is  specified,  as  well 
as  appropriate  values  of  PINF  and  TINF. 


REAL  GAS  ( IRG  =  1) 


All  quantities  in  real  gas  calculations  are  dimensional  with  the 

units: 


p  1 b  f / ft 2 

p  slugs/ft3 

T  °R 


(u.v.w) 


ft/sec 


The  real  gas  thermodynamic  properties  used  in  3DSAP  are  for  equilibrium 
air.  The  gas  properties  required  during  the  calculation  are  obtained  by 
interpolation  on  the  thermodynamic  property  tables  supplied  with  the  code. 
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C.  OUTPUT  CONTROLS 


This  section  details  the  input  variables  that  control  the  type 
of  output  generated  by  the  30SAP  code. 

1.  STANDARD  FORMAT 


The  standard  output  of  3DSAP  consists  of  the  following: 

a. )  surface  pressure  and  shock  shape  for  each  surface  on 

the  nose, 

b. )  summary  of  force  and  moment  data  on  the  nose  at  each 

axial  station  specified, 

c. )  surface  pressure  and  shock  shape  for  each  plane  on 

the  frustum, 

d. )  summary  of  force  and  moment  data  on  the  frustum  at 

each  axial  station  computed. 

Force  and  moment  coefficients  are  normalized  using  a  local  reference 
area  (Trr2)  ard  the  local  reference  length  (z,  measured  from  the  arbitrary 
point  z  =  0).  Options  are  available  to  the  user  to  specify  a  reference  area 
(AREF)  and  a  reference  length  (ALREF)  to  be  used  at  all  stations.  In  addi¬ 
tion,  an  arbitrary  moment  reference  point  may  be  used  by  specifying  XBAR, 

YBAR,  and/or  ZBAR. 

2.  OPTIONAL  OUTPUT 


Additional  output  which  may  be  generated  includes: 

a. )  tables  of  shock  properties  (KENSE(l)  t  0), 

b. )  field  data  (data  printed  every  I  PRINT  stations 

beginning  with  station  MPRINT). 
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D.  SPECIAL  PROVISIONS  AND  OTHER  CONTROLS 


3DSAP  has  input  parameters  and  options  other  than  those  described 
above,  designed  to  provide  the  user  with  a  high  degree  of  flexibility  in 
the  application  of  the  code.  These  parameters  and  options  are  described 
below. 


Global  Transonic  Iteration 


As  discussed  in  References  1  and  2,  3DSAP  has  the  option  of  per¬ 
forming  a  global  iteration  on  shock  slopes  in  the  approximate  transonic 
solution  on  the  nose.  Although  not  normally  exercised  in  the  standard 
operating  procedure  of  3DSAP,  this  global  iteration  can  be  implemented  by 
setting  ISHK  to  zero.  The  use  of  this  global  iteration  is  discussed  in 
Reference  2. 

Surface  Entropy  Relaxation 

In  the  standard  operating  procedure  of  3DSAP  the  surface  entropy 
is  assigned  the  stagnation  value,  since  the  stagnation  streamline  wets  the 
body  surface  in  inviscid  flow.  However,  as  the  length  of  the  frustum 
approaches  infinity,  the  vortical  layer  induced  by  a  blunt  nosetip  becomes 
singular.  (Such  a  situation  would  be  encountered  if  a  sharp  cone  solution 
were  sought  as  the  limit  of  a  blunted  cone  solution  as  the  bluntness  ratio 
goes  to  zero.)  To  avoid  this  problem,  the  surface  entropy  in  the  windward 
plane  of  the  frustum  may  be  relaxed  by  determining  its  value  by  extra¬ 
polation  of  the  entropies  at  the  two  adjacent  field  points  on  the  body 
normal.  This  option  is  exercised  by  setting  ICS  equal  to  zero. 

Termination  of  the  Calculation 


In  normal  usage,  3DSAP  terminates  when  the  end  of  the  frustum 
has  been  reached  in  the  forward-marching  computational  procedure.  For  code 
check-out  purposes,  it  is  sometimes  convenient  to  compute  only  a  predetermined 
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number  of  steps  on  the  frustum.  3DSAP  can  be  terminated  before  the  end  of 
the  frustum  is  reached  by  specifying  the  number  of  computational  steps  to 
be  taken  on  the  frustum,  MTERM. 

Intervals  Across  the  Shock  Layer 

Since  it  is  a  complete  flow  field  technique,  3DSAP  determines 
properties  across  the  shock  layer  (along  body  normals)  as  well  as  surface 
pressures  and  shock  shape.  Two  input  parameters  are  available  in  3DSAP 
to  control  the  point  spacing  between  the  body  and  the  shock  along  body 
normals,  NMAX1  and  NMAX2.  NMAX1  is  the  number  of  points  used  across  the 
shock  layer  in  the  transonic  solution  and  is  usually  set  to  five.  Points 
in  the  transonic  region  are  equally  spaced  in  stream  function  across  the 
body  normal  (vonMises  coordinates).  For  supersonic  calculations,  both  on 
the  nose  and  the  frustum,  the  number  of  points  used  is  determined  by  NMAX2. 
Points  in  the  supersonic  calculation  are  equally  spaced  in  distance  from 
the  body,  measured  along  the  body  normal.  The  standard  value  of  NMAX2  is 
ten,  although  other  values  may  be  appropriate  depending  on  the  freestream 
Mach  number,  as  discussed  in  Section  6. 

Definition  of  Transonic  Region 

3DSAP  automatically  determines  the  end  point  of  the  transonic 
region  in  each  nose  surface  being  computed  by  taking  the  more  downstream 
of: 

a. )  the  body  point  that  is  NSON  input  points  downstream 

of  the  sonic  point  (which  is  determined  from  a  correla¬ 
tion  based  on  local  body  slope),  or 

b. )  the  most  forward  point  such  that  no  downstream  point 

has  a  surface  slope  relative  to  the  freestream  velocity 
vector  of  greater  than  45°. 

For  most  applications,  a  value  of  NSON  equal  to  one  will  provide  an  adequate 
definition  of  the  transonic  regions.  By  careful  choice  of  NSON.  however,  the 
user  has  the  ability  to  enlarge  the  transonic  region,  if  so  desired. 
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Noset i p  Pressure  Correlation  Option 


As  described  in  Section  2  of  this  report,  a  new  nosetip  pressure 
correlation  which  was  developed  for  indented  nosetip  geometries  has  been 
added  as  an  option  to  3DSAP.  This  new  correlation  will  automatically  be 
used  for  arbitrary  nosetip  geometries,  or  if  the  user  sets  IPRESS  =  2. 

The  original  correlation  described  in  Reference  1  will  automatically  be 
used  for  analytic  nosetip  geometries,  or  if  the  user  sets  IPRESS  =  1. 

Normal  Coordinate  Stretching 

The  stretching  of  the  normal  coordinate  to  concentrate  mesh  points 
near  the  vehicle  surface,  as  described  in  section  3,  is  invoked  by  specifying 
a  value  for  the  stretching  parameter,  a*(ALPHSTR).  The  effect  of  this 
stretching  parameter  is  illustrated  in  Figure  3.1.  Different  values  of  a* 
may  be  assigned  to  the  bow  and  embedded  shocks  by  specifying  ALPHSTR(l) 
and  ALPHSTR(2) ,  respectively. 

Pressure  Profile  Exponent 

The  assumed  pressure  profile  across  the  shock  layer  in  the  tran¬ 
sonic  solution  of  3DSAP  takes  the  form 

oM  -  Pb  Mi)ANEXP(Ps-Pb> 

where  ANEXP  =  2  in  the  original  version  of  3DSAP  (Reference  1).  In  the 
current  version  of  3DSAP,  ANEXP  has  been  pre-set  to  0.75  (recommended  value), 
but  the  user  c;;.  modify  this  value  by  specifying  a  different  value  for  ANEXP. 
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APPENDIX  B.  GLOSSARY  OF  INPUT  VARIABLES 


(Default  values  shown  In  parentheses  at  right) 


AB 

ALPHA 

ALPHSTR(I) 


ALREF 

ALT 

AMINF 

ANEXP 

AREF 

BETA 

DELZ 

DZF(I) 


Ratio  of  major  to  minor  axes  for  elliptic 
nosetip  (analytic  elliptic  nosetip  geometry 
option) 

a,  Angle  of  attack  (degrees) 

a*,  Stretching  parameters  for  normal  grid  (0.) 

stretching;  I  =  1  for  Region  I  and  I  =  2 
for  Region  II  (for  flared  bodies)  of  frustum 
calculation 

Reference  length  for  moment  coefficients  (if 
zero,  local  z  used) 

Altitude  (feet),  if  IATMP  used 

Freestream  Mach  number  (if  QINF  not  input) 

Exponent  for  assumed  pressure  profile  across  (0.75) 

the  nosetip  shock  layer  (transonic  region) 

Reference  area  for  force  and  moment  coefficients 
(if  zero,  local  mr^,2  used) 

3,  Sideslip  angle  (degrees) 

Axial  location  of  the  geometric  stagnation  (0.0) 

point  for  a  spherical  nosetip 

Array  of  axial  lengths  of  conical  frustum 
segments  (NCON  values) 


Dl( K)  dj.  Vertical  offset  of  nose  cross-section  center, 

maximum  value  of  K  =  25 


D2(K)  d2,  Lateral  offset  of  nose  cross-section  center, 

maximum  value  of  K  =  25 

GIDEAL  Y.  Ideal  gas  isentropic  exponent  (1.4) 

IATMP  Atmospheric  freestream  property  assignment  (0) 

indicator: 

=  62  Freestream  pressure  and  temperature 
from  1962  standard  atmosphere 

=  59  Freestream  pressure  and  temperature 
from  1959  ARDC  tables 
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IB  IC 

ICS 

I  ELL 
IPRESS 

IPRINT 

IRG 

ISA 

ISHK 

KENSE(l) 

KMAX 

MPRINT 

MSYM 

MTERM 

NCON 

NITER 


f  0  For  analytic  biconic  nosetip  geometry 
option 

Surface  entropy  relaxation  indicator:  (1) 

■  0  Surface  entropy  relaxed  in  wind  plane 
=  1  No  entropy  relaxation 

f  0  For  analytic  elliptic  nosetip  geometry 
option 

Nosetip  pressure  correlation  option  flag 
a  l  For  spherical  nosetip  correlation 
=  2  For  indented  nosetip  correlation 

Print  increment  control  for  field  output  (1) 

Real  gas  option  indicator: 

=  0  For  ideal  gas  thermodynamics 
=  1  For  equilibrium  air  thermodynamics 

Plane  of  symmetry  indicator:  (1) 

=  1  Pitch  plane  of  symmetry  and  S  =  0 
=  2  Asymmetric  shape  or  B  t  0 

Global  transonic  iteration  indicator:  (1) 

=  0  For  global  iteration  on  transonic  shock  slopes 
=  1  For  assignment  of  transonic  shock  slopes 

Shock  table  output  indicator:  (0) 

>  0  Shock  tables  printed  out  as  part  of  calculation 
<  0  Code  generates  shock  tables  only 

Number  of  points  per  meridional  plane  on  nose  -  (15) 

maximum  of  25  (required  for  analytically  defined 
spherical  nose  only) 

Station  at  which  printing  of  field  output  begins  (999) 

Index  of  meridional  surface  about  which  the 
geometry  is  symmetric 

Number  of  steps  to  be  taken  (if  termination  (999) 

desired  before  z  a  ZF(NC0N  +1)) 

Number  of  conical  segments  on  the  frustum  (1) 

(maximum  of  5) 

Maximum  number  of  global  transonic  iterations,  (100) 

if  ISHK  =  0 


I 


NMAX1 

NMAX2 

NSON 

NUWTB 

PHI(M) 

PHIBD(M) 

PINF 

QINF 

RBD( K,M) 

RC 

RF(  I) 

R I DEAL 
RN 

STAB 
TEST ( 11 ) 
TEST( 12) 
TEST(28) 


Points  between  body  and  shock  -  transonic  (5) 

(maximum  of  20) 

Points  between  body  and  shock  -  supersonic  (10) 

(maximum  of  20) 

Number  of  points  beyond  sonic  point  to  be  (0) 

included  in  transonic  region 

Number  of  entries  in  shock  table  (91) 


4>,  Meridional  angles  of  surfaces  to  be  ( ISA  =«  1  0°, 90°.  180°) 
calculated  on  nose  (deg.)  ( ISA  =  2  0°, 90°, 180°, 270°) 

(maximum  value  of  M  =  8) 


Meridional  angles  of  surfaces  in  which 
arbitrary  nose  geometry  is  defined  (Deg.) 

(maximum  value  of  M  =  8) 

p^,  Freestream  static  pressure  (Real  Gas,  1 bf/ft2 )  (1.0) 

qoo,  Freestream  velocity  (if  AMINF  not  input) 

(Real  Gas,  ft/ sec) 

r,  Local  nose  radius  at  axial  station  K  and 
meridional  surface  M;  maximum  value  of  K  *  25, 
maximum  value  of  M  =  8 

Comer  radius  for  biconic  nosetip  geometry 

Array  of  radii  of  conical  frustum  segments 
(NC0N  +  1  values) 


R,  Gas  constant  (Ideal  Gas)  (1.0) 

r^j,  Nose  radius  (spherical  or  biconic  noses) 

Stability  factor  on  CFL  condition  for  super-  (1.0) 

sonic  calculations 

Convergence  criterion  for  shock  property  (10-“) 

calculations  (Real  Gas) 

Convergence  criterion  for  stagnation  condition  (10'“) 

iteration  (Real  Gas) 

Test  for  initiation  of  bow  and  embedded  shock  (0.05) 

intersection  procedure;  shock  intersection 
assumed  when  (ys^-ys2)/(ys^)  *  TEST(28) 


» 
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TEST( 29) 
TEST( 30) 

TEST (31) 

TEST ( 32 ) 

TEST (33) 

TEST ( 34 ) 

TEST( 35) 

TEST(36) 

TEST(37) 

TEST( 38) 

TEST (39) 

TEST(40) 

THF(I) 

THN 

TINF 

XBAR.YBAR, 

ZBAR 

ZBD(K) 

ZF(  I) 

Z I DEAL 


Convergence  criterion  for  embedded  and  bow  shock  (0.001) 
intersection  Iterative  calculation 

Test  for  initiation  of  flare  calculation;  flare  (2.0) 
assumed  when  THF( I  +  1)  -  THF(I)  >  TEST(3Q), 
in  degrees 

Convergence  criterion  for  global  transonic  (0.005) 

iteration 

Damping  factor  for  global  transonic  (0.6) 

iteration 

Convergence  criterion  for  local  transonic  (0.005) 

iteration 

Convergence  criterion  for  supersonic  shock  (10“5) 

calculation 

Alternate  convergence  criterion  for  supersonic  (10~“) 

shock  calculation 

Acceleration  factor  for  local  transonic  (1.1) 

iteration 

Convergence  criterion  for  flare  embedded  shock  (10“4) 

initialization  iteration 

Increment  in  axial  location  used  in  Initializing  (0.05) 
embedded  shock  calculation 

Criterion  for  respacing  grid  points  in  region  (0.05) 

between  embedded  and  bow  shocks 

Criterion  for  respacing  grid  points  in  region  (0.05) 

between  body  and  embedded  shock 


Array  of  cone  angles  of  conical  frustum 
segments  (degrees)  (NC0N  values) 

Nose  cone  angle  for  biconic  nosetlps 

T^,  Freestream  static  temperature  (Real  Gas,  OR)  (1.0) 

_  __ _  i 

(X,Y,Z)  Moment  reference  point 


Axial  stations  for  specification  of  nose 
geometry,  maximum  value  of  K  s  25 

Array  of  axial  stations  of  conical  frustum 
segments  (NC0N  +  1  values) 

z.  Compressibility  factor  (Ideal  Gas)  (1.0) 


49. 


DISTRIBUTION  LIST 


Ballistic  Missile  Office 
BMO/SYOT 

Attn:  Maj.  K.  Yelmgren  (2) 

Norton  AFB,  CA  92409 

Defense  Technical  Information  Center  (2) 
Cameron  Station 
Alexandria,  VA  22314 

Air  University  Library 
Maxwell  AFB,  AL  36112 

TRW  OS SO 

Attn:  W.  Grabowsky  (2) 

P.  0.  Box  1310 

San  Bernardino,  CA  92402 

TRW  Systems  Group  (2) 

Attn:  0.  Ohrenberger 
M.  Gyetvay 
1  Space  Park 
Redondo  Beach,  CA  92078 

Headquarters,  Arnold  Engineering 
Development  Center 
Arnold  Air  Force  Station 
Attn:  Library/Documents 
Tullahoma,  TN  37389 

Armament  Development  and  Test  Center 
Attn:  Technical  Library,  DLOSL 
Eglin  AFB,  FL  32542 

Air  Force  Wright  Aeronautical  Laboratories 
Air  Force  Systems  Command 
Attn:  M.  Buck  (AFWAL/FIM) 

R.  Neumann  (AFWAL/FIMG) 

V.  Dahlem  (AFWAL/FIMG) 
Wright-Patterson  AFB,  OH  45433 

U.  S.  Army  Ballistic  Missile 
Defense  Agency/ ATC-M 
Attn:  J.  Papadopoulos 
P.  0.  Box  V500 
Huntsville,  AL  35807 

Director,  Defense  Nuclear  Agency 
Attn:  J.  Somers  (SPAS) 

Washington,  DC  20305 


Naval  Surfact  Weapons  Center 
Attn:  Carson  Lyons/ KQ6 
White  Oak  Laboratories 
Silver  Spring,  MD  20910 

Acurex  Aerotherm 
Aerospace  Systems  Division 
Attn:  C.  Nardo 
485  Clyde  Avenue 
Mountain  View,  CA  94042 

Avco  Systems  Division 
Attn:  N.  Thyson 
201  Lowell  Street 
Wilmington,  Mass  01887 

General  Electric  Company 
Attn:  R.  Neff 
3198  Chestnut  Street 
Philadelphia,  PA  19101 

Lockheed  Missiles  and  Space  Co. 

P.  0.  Box  504 

Attn:  G.  T.  Chrusciel 

Sunnyvale,  CA  94086 

McDonnell  Douglas  Astronautics  Co. 
Attn:  J.  Copper 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

PDA  Engineering 
Attn:  M.  Sherman 
(3)  1560  Brookhollow  Drive 
Santa  Ana,  CA  92705 

Sandia  Laboratories 
P.  0.  Box  5800 
Attn:  Library 
Albuquerque,  NM  87115 

Science  Appl ications,  Inc. 

Attn:  A.  Martel lucci 
994  Old  Eagle  School  Road 
Suite  1018 
Wayne,  PA  19087 


